Pyramid diffraction in parity-time-symmetric optical lattices 
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Nonlinear dynamics of wave packets in two-dimensional parity-time-symmetric optical lattices 
near the phase-transition point are analytically studied. A novel fourth-order equation is derived 
for the envelope of these wave packets. A pyramid diffraction pattern is demonstrated in both the 
linear and nonlinear regimes. Blow-up is also possible in the nonlinear regime for both focusing and 
defocusing nonlinearities. 



Parity-time ('PT)-symmetric wave systems have the 
unintuitive property that their linear spectrum can be 
completely real even though they contain gain and loss 
PQ. In spatial optics, PT-symmetric systems can be re- 
alized by employing symmetric index guiding and an an- 
tisymmetric gain/loss profile [2-5 . In temporal optics, 
T^T-symmetric systems can be obtained as well [SHE]- So 
far, a number of novel phenomena in optical VT sys- 
tems have been reported, including phase transition, non- 
reciprocal Bloch oscillation, unidirectional propagation, 
distinct pattern of diffraction, formation of solitons and 
breathers, wave blowup, and so on [4HT4]. 

In this Letter, we analytically study nonlinear dynam- 
ics of wave packets in two-dimensional PT-symmetric 
optical lattices near the phase-transition point (where 
diffraction surfaces of Bloch bands cross like the inter- 
section of four planes). Near these intersections we show 
that the evolution of wave packets is governed by a novel 
fourth-order equation. Based on this envelope equation, 
we predict a pyramid (i.e., expanding square) diffraction 
pattern in both linear and nonlinear regimes. Further, in 
the nonlinear regime blow-up can occur for both focus- 
ing and defocusing nonlinearities. These predictions are 
verified in the full equation as well. 

The model for nonlinear propagation of light beams in 
PT-symmetric optical lattices is taken as 



Diffraction Surfaces 



Envelope Solution 



i* z + V 2 * + V(x, y)* + ct|*| 2 * = 0, 



(1) 



where z is the propagation direction, (x, y) is the trans- 
verse plane, V 2 = $ 2 + <9 2 , and a = ±1 is the sign of non- 
linearity. The T^T-symmetric potential V(x, y) is taken 
as V{x, y) = V(x) + V(y), where 



V(x) = V 2 [cos(2x) + iW sm(2x)] , 



(2) 



Vq is the potential depth and Wo is the relative gain/loss 
strength. 

We begin by considering the linear diffraction rela- 
tion of Eq. at the phase-transition point Wo = 1 
[9l [13]. In this case, the linear equation of ([I]) can 
be solved exactly |T3 . The diffraction relation is fi = 
(k x + 2mi) 2 + (k y + 2m2) 2 , where (k x ,k y ) are Bloch 
wavenumbers in the first Brillouin zone 1 ^ k xi ky ^ 1, 
and (mi, 777,2) are any pair of nonnegative integers. The 
most complex degeneracies occur at points k x = 0, ±1 
and k y = 0,±1, where the diffraction surface intersects 
itself four-fold as illustrated in Fig. [I] If a carrier Bloch 




FIG. 1: (Upper left) Diffraction relation near the intersection 
point (k x ,k yi fi) = (1,1,2) (marked by a red dot). (Upper 
right) Linear-diffraction pattern of an initial Gaussian enve- 
lope at phase transition in the envelope equation ( 11 ). (Lower 
row) Linear diffraction of an initial Gaussian wave packet at 
phase transition in the full equation |l]). 



wave is chosen at one of these degeneracies, then the enve- 
lope of the resulting wave packet exhibits novel behavior 
which we elucidate below. 

To analyze the structure of these degenerices we exploit 
the fact that the potential V(x, y) is separable. Posed as 
an eigenvalue problem for \I/ = cj){x, y)e~ lflz in the linear 
equation of |l]), we get L<j) = — /i</>, where L = -\-L^ y \ 
= <9 2 + V (x), and V (x) is the VT lattice Q at the 
phase-transition point Wo = 1. At four- fold intersection 
points, the eigenvalues are /i = n\ + 772, where (ni, 772) 
are any pair of positive integers. The operator (L^) 
has eigenvalues n\ (77?,) with geometric multiplicity 1 and 
algebraic multiplicity 2 [14]. Let <j) ei (x) {4> e<2 {y)) be the 
eigenfunction and (j) 9l {x) ((j) g2 {y)) the associated gener- 
alized eigenfunction. Then 



^(x) = I ni (Voe ix ), r 2 (v) = LAVoe'y), (3) 

where I n (yoe lx ) is the modified Bessel function J n (Vbe 2:c ) 
normalized to have unit peak amplitude, and (L^ + 

n 2 )^ 1 = (j) e \ (L^+n 2 ,)^ 2 = cj) e \ Since L = L^+L^\ 
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FIG. 2: Pyramid diffraction of a Gaussian wavepacket in the 
linear equation below phase transition. Upper row: diffrac- 
tion in the envelope equation (11); lower row: diffraction in 
the full equation |l]). 



we see that L has two eigenfunctions 

4P 1 {x,y) = r i (xW 2 (y), (4a) 
<j> 02 (x,y) = 4> ei (x)<j> 92 (y) ~ 4> 9l (x)4> e2 (y). (4b) 

In addition, the first eigenfunction </> 01 has two general- 
ized eigenfunctions 

^Hx : y) = [r i (x)^(y)^^(x)r 2 (y)}/2, (5a) 
4> 21 {x,y) = ^{x)^{y), (5b) 

where (L + /i)^ 11 = ^ 01 , and (L + /i)^ 21 = 

We now study the nonlinear dynamics of wave packets 
near these intersections. For simplicity, we will conduct 
the analysis at the lowest intersection point, \i = 2. But 
similar results can be obtained for any intersection point. 

The perturbation expansion for the wave packet near 
the intersection \i = 2 is 



\I> = e 2 ^e 



-i/it 



ip = i/j + e0i + . . . , 



where iJjq = A(X,Y, Z)(j) 01 (x,y) is the leading-order 
wave packet for the Bloch mode </> 01 at the intersec- 
tion point (k x ,k y ,fi) = (1,1,2), (X,Y,Z) = (ex : ey,ez) 
are slow spatial variables, and < e <C 1. Near the 
phase-transition point, Wo can be expressed as Wo = 
1 — r]€ 2 /Vq, where r] measures the deviation from phase- 
transition. After introducing the slow variables into 
equation ([!]), the new equation for ip is 



(L + /i)V> = -ie^z - 2e(il> xX + ^ yY ) - e 2 (d 2 x H 
+ irj€ 2 [sm(2x) + sin(2j/)]^ - e 3 cr|^| V- 



(7) 



We proceed by inserting expansion (|6| into equation 
^ and solving for ^ n at each order. Each ijj n satisfies 
a linear inhomogeneous equation, with the homogeneous 
operator being L + /i. In order for it to be solvable, 
the Fredholm conditions need to be satisfied, i.e., the 



inhomogeneous term must be orthogonal to the kernels 
and (/> 02 * of the adjoint operator L* + fi. Here * 
stands for complex conjugation. 

At 0(e) the solvability conditions for ipi are automat- 
ically satisfied, and thus we can solve ipi as 



-I An 



2A> 



2 A 



±02 



(8) 



where B(X, Y, Z) is the envelope of the second eigenfunc- 
tion ^ 02 , (L + ii)(\) a = 0°\ (L + /i)/ = 0°\ and a , b 
are assumed to be orthogonal to </> 01 and 02 . 

Now we proceed to the ^2 equation at 0(e 2 ). The or- 
thogonality condition with is automatically satisfied, 
and the orthogonality condition with 02 * gives 



iB z 



2 A 



XX 



2A 



YY, 



(9) 



which defines the connection between envelopes A and B 
of the two eigenmodes at the Bloch-surface intersection. 
Under this relation, the ^2 equation can be solved. 

Finally we proceed to the ^3 equation at 0(e 3 ). The 
orthogonality condition with 



i01* 



gives 



d%A 



8(d 2 



x 



d\)d z A 



+8{d 2 x -d 2 Y )(iB) 



■ x -d*){d x - Oy)a 

ad z A + ia\A\ 2 A = 0, (10) 



where 



W \ 



^ \^\ 2 4> 01 4> 21 dxdy 



Jo Jo ' 



,01^21 



dxdy 



This equation, combined with equation (|9|, yields a sin- 
gle fourth-order envelope equation for A as 



8(d 2 x + d\)d\A + 16(c>| 
ad 2 z A + \dd z {\A\ 2 A) = 0. 



d 2 fA 



(11) 



This novel envelope equation is one of the main results 
in this Letter. 

It remains to relate the initial conditions for ^ with 



(6) those for the envelope equation (11). By collecting the 



^0,^1 an d ^2 solutions from the above analysis and 
projecting the resulting perturbation-series solution ([6| 
onto the eigenfunctions and generalized eigenfunctions, 
we find the dominant terms of \I/ are given by 



01 



/02 



where 



D 



C = -i(A 2 



2 A 



x 



■2A Y ), 



(12) 



(13) 



- A zz - 2A ZX - 2A ZY - -A + 4iB x - 4LB y . (14) 



Thus, from initial envelope functions A,B,C,D of the 
eigenfunctions and generalized eigenfunctions in we 
can obtain initial conditions for A, A Zl Azz and Azzz 



from (10), (13) and (14). 
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FIG. 3: Nonlinear dynamics of wave packets below phase 
transition. Upper row: env elope solutions in (11) at Z « 2 
for three values of Ao in ( 15 ). Lower row: solutions of the full 
equation |l]) for the initial wavepacket with Ao = 6 (left) at 
later distances under focusing (middle) and defocusing (right) 
nonlinearities. 



Direct simulations show strong agreement between the 
envelope dynamics in (11) and those corresponding wave 
packets in the full equation ([!]). For this letter we take 
the initial conditions 



A e 



~{X 2 +Y 2 ) 



0, & 



%A 



-ia\A\ 2 A 
(15) 

in the envelope equation, or the equivalent initial condi- 
tions based on (12) for simulations of the full equation 



([!]). For the constants we take V§ = 6, e = 0.1, and 77 = 
or 1 (at or below phase transition respectively). Then we 
find a = 12^, and a « 7.3a. 

In the linear equation at the phase transition point, 
i.e., a = a = 0, equation ( pTj ) has the general solution 

A = A 1 (X - 2Z,Y - 2Z) + A 2 (X — 2Z, Y + 2Z) 
+ A 3 (X + 2Z, Y - 2Z) + A 4 (X + 2Z, Y + 2Z) (16) 

for arbitrary A n functions. In general this corresponds to 



an expanding square wave front propagating with speeds 
±2 in both X and Y directions, which we term pyramid 
diffraction. This pattern is illustrated in Fig. [I] for both 
the envelope and full equations. Notice that the wave 
fronts are flat on all four sides. 

In the linear equation but below the phase transition 
point (a = 12, a = 0), the pyramid diffraction is qualita- 
tively similar to that in Fig. 1, with wave fronts expand- 
ing roughly like a square. But the wave fronts are no 
longer flat. In addition, the core develops an 'x' shape. 
An example is shown in Fig. [2j where diffractions in both 
the envelope and full equations are displayed. 

In the presence of nonlinearity (a « 7.3a) and be- 
low phase transition, the wave packet diffracts away if 
its initial amplitude is below a certain threshold value. 
This nonlinear diffraction is also pyramid-like, closely re- 
sembling the linear pyramid diffraction in Fig. 2. An 
example is displayed in Fig. 3 (upper left panel). How- 
ever, if the initial amplitude is above this threshold, the 
envelope solution blows up to infinity in finite distance. 
For example, with the initial condition (15), the envelope 
solution in (11) blows up when Aq > 3.2. These blowup 



solutions are displayed in Fig. 3 (upper middle and right 
panels). Remarkably, this blowup is independent of the 
sign of the nonlinearity, a fact which is clear from the 
envelope equation (11), since a sign change in a can be 



accounted for by taking the complex conjugate of this 
equation. In the full equation 0, we have confirmed 
that similar growth occurs for both signs of the nonlin- 
earity as well. For instance, evolution of a wavepacket 
(corresponding to Ao = 6) under focusing and defocus- 
ing nonlinearities are displayed in Fig. 3 (lower row). In 
both cases solutions rise to very high amplitudes as the 
envelope equation predicts. 

In summary, we have analyzed nonlinear dynamics of 
wave packets in two-dimensional T^T-symmetric lattices 
near the phase-transition point. In the linear regime, 
pyramid diffraction is demonstrated. In the nonlinear 
regime, wave blowup is obtained for both focusing and 
defocusing nonlinearities. 
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